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A quantum linear Boltzmann equation is proposed, constructed in terms of the operator-valued 
dynamic structure factor of the macroscopic system the test particle is interacting with. Due to this 
operator structure it is a non-Abelian linear Boltzmann equation and when expressed through the 
Wigner function it allows for a direct comparison with the classical one. Considering a Brownian 
particle the corresponding Fokker-Planck equation is obtained in a most direct way taking the limit 
of small energy and momentum transfer. A typically quantum correction to the Kramers equation 
O^l \ thus appears, describing diffusion in position and further implying a correction to Einstein's diffusion 

coefficient in the high temperature and friction limit in which Smoluchowski equation emerges. 
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Oh The study of irreversible processes, even in relatively simple cases in which one considers the reduced dynamics of a 
| system with few degrees of freedom interacting with a suitable reservoir, is a subject of major interest in both classical 
and quantum mechanics and a necessity whenever one wants to cope with real experiments, possibly starting from 
the underlying microphysical dynamics. In particular in quantum mechanics irreversible incoherent, and possibly 
dissipative, processes play a decisive role in determining the feasibility of quantum computers [Q and under the 
i-H ' unifying concept of decoherence might even have a deep impact in the understanding of quantum mechanics and of its 
\ relationship with the classical description of the world [^) , provided a precise connection with microphysical models 
\l ■ can be given. 

In the specific case of a microsystem whose interaction with a macroscopic system or reservoir can be described 
in terms of collisions, a major tool for the classical description of its dynamics is the linear Boltzmann equation 
or transport equation. This equation has proved essential in fields like neutron transport theory |3) and is still the 
, standard reference for this kind of dynamics even at quantum level ^ , the only improvement lying in the introduction 
of the quantum collision cross-section in place of the classical one. This kinetic equation, provided a suitable expansion 
of the collision term can be performed, then leads to a Fokker-Planck equation for the study of dissipation H , which 
for Brownian motion is the so-called Kramers equation. Considering situations in which dissipation is very important, 
Kramers equation in turn leads to the so-called Smoluchowski equation in which Einstein's diffusion coefficient appears. 

In this paper we give the structure of a fully quantum, non-Abelian linear Boltzmann equation, non-Abelian just due 
to the presence in it of the operator- valued rather than c-number dynamic structure factor of the macroscopic system, 
as typical of the quantum realm, giving a physical example of a recently obtained mathematical result on Lindblad 
generators of translation-covariant completely positive quantum dynamical semigroups || . Considering the Brownian 
, motion of a test particle in a gas a most straightforward expansion of the non-Abelian linear Boltzmann equation, 
written in terms of the Wigner function, leads to a Kramers equation with a typical quantum correction linked to 
r> \ position diffusion vanishing in the semiclassical limit fi — > 0. Such a correction ensures positivity at quantum level and 
is connected to the time evolution of the off-diagonal matrix elements of the statistical operator in the momentum 
representation, which can only be given in a fully quantum mechanical version of the linear Boltzmann equation. 
Studying the high friction limit of this Fokker-Planck equation analogously to the classical case, in the same spirit, 
also a quantum correction to Einstein's diffusion coefficient appears. 

I. NON-ABELIAN LINEAR BOLTZMANN EQUATION 

Let us start writing the linear Boltzmann equation for the classical case, still relying on the quantum expression 
for the collision cross-section, in a way which is actually not the most common one. According to a famous paper by 
van Hove , the energy-dependent differential cross-section per target particle describing scattering of a microscopic 
probe off a macroscopic sample is given by 

P -Z(q)S(q,E) (1) 



dVLpi dEpi p 

if the momentum of the microscopic probe changes from p to p' = p + q. is here the collision cross-section 

for the single scattering event, given by = M 2 (2irTi) 3 ^-\t(q)\ 2 , with i(q) Fourier transform of the T-matrix 
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describing the one-pair collisions between the test particle and the particles constituting the sample, supposed to 
depend only on the modulus of the momentum transfer q. S(q, E) is a two-point correlation function known as 
dynamic structure factor H, reflecting the equilibrium many-body properties of the fluid, depending on momentum q 
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and energy E = E(q, p) = j^j + ^f- transferred to the test particle of mass M, whose general expression is given by 
the Fourier transform with respect to energy and momentum transfer of the time dependent spatial autocorrelation 
function according to: 



S(q,E) = ^LJdtJ d 3 xe*( Et - q -^ 

x^Jd 3 y (N(y)N(x + y,t)). 



Exploiting (Q) the linear Boltzmann equation in the absence of external potentials may be written as Q 

9 f( x ,p,t)=-^-.V x f(x,p,t) + ^ I 'd 3 qX(q) (2) 



dt Jy ' M '^v-^v M2 

x [S(q,p- q)f(x,p- q,t) - S(q,p)f{x,p,t)} 

with S(q,p) = S(q,E) the dynamic structure factor of the homogeneous macroscopic fluid of density n the particle 
is interacting with. So far the classical case, all expressions appearing in (^|) being c-number functions. 

Recent work on the study of subdynamics in non-relativistic quantum field theory |1 has led to a scattering theory 
derivation of a particular structure of master-equation JIo|-^2[, which goes beyond the usual limitation of linear 
coupling and is a natural candidate to be considered as a fully quantum linear Boltzmann equation, as we will show in 
the sequel. Instead of an equation for the classical distribution function it is an equation for the full statistical operator 
p associated to the test particle, describing not only its momentum distribution, but also its coherence properties. 
These are particularly relevant for the time evolution of the off-diagonal matrix elements, related to the phenomenon 
of decoherence, and for the low temperature behavior, where the classical particle viewpoint is insufficient and wave 
effects are to be taken into account. The quantum kinetic equation is given by: 

I — iWHfl + ^/AlX.) (3) 



with x and p position and momentum operator for the test particle, Ho = p 2 /2M, S(q, p) the operator-valued 
dynamic structure factor, which is a positive operator due to (Q). Equation ([|) addresses the same physics as 
the linear Boltzmann equation but is non-Abelian due to the appearance in it of operators instead of c-number 
functions and it provides a physical realization of the general mathematical structure of translation-covariant quantum 
dynamical semigroup other known physical examples restricting to the diffusive case. Neglecting the dependence 
on the momentum operator p and therefore on the energy transfer of the dynamic structure factor in (||), so that 
its appearance only corrects the integration measure, one recovers a structure of master-equation typically proposed 
on a more phenomenological basis for the study of the phenomenon of decoherence, which has gone through many 
refinements |l3| , all missing the correct description of momentum and energy transfer, which may not be critical for 
the description of short-time decoherence but are crucial for the approach to the correct stationary state. In fact 
while the models in |l3f| predict a steady growth of kinetic energy and do not admit a stationary state, provided the 
fluid is in a /3-KMS state or equivalently the dynamic structure factor satisfies the detailed balance condition, eq. (||) 
admits a stationary solution of the form exp(— /3p 2 /2M), with M mass of the test particle and inverse temperature 
of the fluid In order to clarify analogies and differences between (||) and (^) we will now consider a dynamics 
in which als the off-diagonal elements of the statistical operator are of relevance, corresponding to position diffusion, 
focusing on a system in which the dynamic structure factor can be explicitly calculated, a free gas of particles of mass 
to satisfying Maxwell-Boltzmann statistics (for the extension to quantum statistics see jl4]]), the expression of the 
correlation function being in this case incidentally the same for both classical and quantum realm [pi 



.V„ l( ^, 1 = A /£_- e -« ? • (4) 
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II. QUANTUM CORRECTION TO KRAMERS EQUATION 



As a first step we recover from (Q) the classical Kramers equation for the Brownian motion of a massive test particle 
in a gas of much lighter ones. Calling a = m/M the ratio between the masses and considering the Brownian limit 
of small a, corresponding to small energy transfer, the correct limiting expression of (^), still satisfying the detailed 
balance condition is |C3 



S ™ (n p) = P-^le-^e-^v) 
ub1?,W V 2tt q 



(5) 



Starting from the linear Boltzmann equation written in terms of the dynamic structure factor as in (g) and exploit- 
ing (||) one obtains in a most straightforward way Kramers equation, a point usually neglected in textbooks ||. 
Substituting (||) in (||) one in fact has 



^f(x,p,t) = -jj ■ Vxf(x,p,t) 



(0) 



t) 



and considering the limit of small momentum transfer q, expanding the exponentials and keeping terms up to second 
order one immediately has for an isotropic medium (qf = hq 2 ) 



^f(x,p,t) = --^ • V x f(x,p,t) 



M d , 

Pi+ JWi ]f{x,p,t) 



(7) 



and therefore Kramers equation for the description of Brownian motion 



^f(x,p,t) = - jj ■ V x f(x,p,t)+ri 
with a friction coefficient explicitly given by 



M 



Vp-(pf(x,p,t)) + —A p f(x,p,t) 



(8) 



We now want to obtain the corresponding result from the non-Abelian linear Boltzmann equation (g). To do this 
it is particularly convenient to use the Wigner function, which even if it is not a well-defined probability density, well 
serves the purpose of comparing quantum versus classical equations. Introducing the Wigner function corresponding 
to the statistical operator p by 



d 3 k 

(2-Khf 



en x ^{p + k/2\p\p-k/2), 



eq. (||) may be easily rewritten for a Brownian particle in a Maxwell-Boltzmann gas, exploiting (^|), as 

d_ 

df 



^/w( a . )Pit ) = _Z.v,./;:i.r.p./, 



M 
n 

'~M 2 V 2vr 



(9) 



e q Vp - cos 



(El 

[am 



e 2M 



The comparison between (^|) and (^) is straightforward: the cosine term in (^) , arising from the non-Abelian structure 
of (^) , is replaced by a factor one in the classical case, and this result can be simply obtained taking the semiclassical 
limit h — » in (0), the argument of the cosine depending linearly on Ti. The quantum correction to Kramers equation 
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then comes from this contribution and amounts to a term corresponding to diffusion in position, as can immediately 
be seen considering the small momentum transfer limit of (^ 

^f?(x,p,t) = -jL ■ V x f?(x,p,t) (10) 



+V 



2 



V P • (pf?(x,p,t)) + jA p f?(x,p,t) + ^A x f?(x, P> t) 



with r\ as in (|g) . Thus position and momentum diffusion do appear together in the quantum description of Brownian 
motion, even though with different weights, as appears from the different (3 dependence of the two contributions, 
the additional diffusion term vanishing in the limit % — > 0. This can be seen more clearly introducing the thermal 

momentum spread Ap 2 th = and the thermal position spread or square thermal wavelength Aa^ h = ^jj, satisfying 

the minimum uncertainty relation Ap th Ax th = so that one has 

p7 = -j] ■ v -/p w + »? l^p ■ W) ( n ) 

The Fokker-Planck equation (|l^) or ([ll]) can be recast in operator form, leading to the well-known structure jl5| 

f -->.a-£!|S>{fc.»] m 

i—l 

7 3 1 3 

~ -jT2 A P^ H fil ~ p4 A ^ h ^ ' 

which could have been obtained directly from (||) as in (lj], without going through the use of the Wigner function, 
missing however in this way the particularly manifest appearance of the quantum correction which provides a inho- 
mogeneous contribution to the collision term in the linear Boltzmann equation through the appearance of the cosine 
in (^|). The appearance of both position and momentum diffusion in ( pi] ) or equivalently ( |l2| ) is due to the fact 
that (§) is an evolution equation not only for the diagonal elements of the statistical operator giving the momentum 
distribution of the test particle, but also for its quantum coherence properties, which actually require preservation 
of the positivity of the statistical operator. In fact the quantum correction in (pi]), actually corresponding to the 
double commutator with the momentum operator in ( |l2|) , is necessary for (|l2| ) to have a Lindblad structure ensuring 
complete positivity (actually corresponding with positivity for Fokker-Planck equations of the form (|l2|) |l6|]). Cor- 
rectness and relevance of this term are heavily debated [p^JlO|]Tl| , p^| , also in connection with the meaning of complete 
positivity |l8| , and its appearance as a typical quantum correction coming from the Brownian and small momentum 
transfer limit of the non-Abelian linear Boltzmann equation (|jj|) expressed in terms of the operator-valued dynamic 
structure factor sheds light on its physical origin. The quantum Fokker-Planck equation (|l^) due to the minimum 
uncertainty relation Ap th Ax th — ^ has furthermore the distinguishing feature of being expressible in terms of a single 
generator for Cartesian direction J]10|, a property actually required for the diffusive component of the generator of a 
translation-covariant semigroup ||7both (|3|) and ( JT^ ) are proper generators of completely positive quantum dynam- 
ical semigroups, admitting the correct equipartition stationary solution and being invariant under translations. The 
three features of complete positivity, equipartition and covariance are in fact not contradictory provided one considers 
covariance under the symmetry relevant to the physical problem, which is not necessarily translation invariance [[l9[ , 
otherwise, independently on positivity or complete positivity of the time evolution, covariance would lead to high 



non-uniqueness of the stationary solution |20[. Equations (|3|) and (12) are essentially the two possible structures 
of generators of translation-covariant semigroups, the results exposed in j(| giving the quantum counterpart of the 
well-known classical result according to which the only meaningful approximation of the integro-differential linear 
Boltzmann equation as a differential equation is a second order Fokker-Planck equation pl[ | . 

Now that the quantum correction to Kramers equation has been given, which unlike previous results is not dependent 
on the presence of external potentials p2[ , but only relates to Heisenberg's uncertainty principle |l5| ], we move on 
to investigate whether this correction also has consequences on the Smoluchowski equation for the description of 
Brownian motion in the limit of high temperature and friction, when momentum is expected to quickly relax to its 
equilibrium value and one is only interested in the time evolution of the slowly varying marginal position distribution 

<rj(x,t) = J d 3 pf™(x,p,t). 
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III. QUANTUM CORRECTION TO SMOLUCHOWSKI EQUATION 



We therefore start from (|iCl|), essentially following the derivation of the Smoluchowski equation from the Kramers 
equation in the high friction limit given by van Kampen |23j] , though more refined derivations have been given pl| , 
differing however only for higher order contributions. Setting 

(13) 



P h *> a 2 _ n 
^16 = 4^ = ^ 



and having in mind a high temperature, strong friction situation in which the quantum position diffusion term can 
be considered as a perturbation, we write (|Io| ) as 



(14) 



— f » _i_ R. . v f w - n a r 

dt p M p LJ xx'-*-x3 p 



Introducing the differential operators 



KM [g] = V p ■ (pg) + jA p g 



-L^ [g] = - 

and expressing as an expansion in powers of 1/rj 



d p 

dt 9+ M ' ^ x9 ~ DxxAx9 



J n J , 



w(0) 



fw(l) 

Jp 



(15) 



eq. ( [LI] ) may be rewritten as [/*] = jjL^ [/*] and solved by iteration using ( [l5"l) and equating on both sides 
contributions of the same order in r\. The only difference with respect to the classical situation considered in p5| ] lies 
in the structure of the iP^ operator, where also the Laplacian with respect to position appears. Going through the 
very same procedure, integrating over the fast degree of freedom given by the momentum dependence of the Wigner 
distribution function, and neglecting terms higher than first order in ~ one obtains again for the marginal position 
distribution a^(x,t) the Smoluchowski equation 



dt a » 



( 1 



\rjM(3 



'(x,t), 



(16) 



with a diffusion coefficient which is however not simply Einstein's, but has a quantum correction, vanishing in the 
semiclassical limit h — > 0, just given by the coefficient (O) responsible for position diffusion in the quantum Kramers 
equation (|Io|). The new overall coefficient can be written as l/(r]M/3)(l + (r][3h) 2 /16), so that the correction is actually 
given by the square ratio between two characteristic times, ph for the bath and 1/rj for the test particle, which should 
be much less than one according to the Markov approximation. 

We have thus proposed a fully quantum, non-Abelian linear Boltzmann equation given by (||) and expressed in 
terms of the T-matrix describing collisions and the operator- valued dynamic structure factor of the environment, thus 
going beyond the usual restriction of linear dissipation. Considering a massive test particle interacting with a gas of 
much lighter ones the linear Boltzmann equation written in terms of the dynamic structure factor leads in a most 
straightforward way to the Kramers equation for the description of Brownian motion in both classical and quantum 
case. One simply has to consider the clear-cut physical limit in which momentum transfer and energy transfer (or 
equivalently ratio between the masses) are small, thus obtaining in the classical case Kramers equation and in the 
quantum case ( |l0| ) with a peculiar, heavily debated quantum correction linked to position diffusion, whose appearance 
is to be traced back to the h dependent, inhomogeneous cosine correction in (|^) vanishing in the semiclassical limit. 
Further studying the case of high temperature and friction the quantum Smoluchowski equation ( |l6| ) is obtained, 
in which Einstein's diffusion coefficient is slightly modified due to the presence of the quantum position diffusion 
coefficient (|l3|), which depends linearly on the inverse temperature and vanishes if h — ► 0. The non-Abelian linear 
Boltzmann equation due to its intrinsic quantum structure should prove as a sound starting point for the study 
of quantum kinetic, in which not only the momentum distribution but also coherence properties are of relevance; 
it furthermore has the advantage of being expressed in terms of quantities of direct physical meaning, such as the 
dynamic structure factor and the collision cross-section, for which suitable phenomenological Ansatz or experimentally 
determined expressions can be inserted. 
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